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Recently, a new class of restricted gravitational wave search templates, termed the Tay- 
lorEt template was proposed for the search of inspiralling compact binaries. The TaylorEt 
approximant is different from the usual time-domain post-Newtonian approximants in 
5-H ' that it employs the orbital binding energy rather than the orbital frequency or the closely 

related parameter "x" . We perform detailed studies to probe the fitting factors of Tay- 
lorEt at 3.5pN for nonspinning comparable mass compact binaries vis-a-vis the TaylorTl, 
TaylorT4, and TaylorF2 at 3.5pN approximants in LIGO, Advanced LIGO and Virgo 
' interferometers. 
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' 1. Templates and LAL routines 

In this article, we consider data analysis performances of several inspiral templates 
for nonspinning compact binaries, used by the LSC Algorithms Library (LAL). 
These routines are compared with the TaylorEt template with respect to the nu- 
merical values of the fitting factor (FF) and the faithfulness [1]. The task is to find 
■ out whether they can be used as effective and faithful search templates if the detec- 

tors strain data is assumed to provide the TaylorEt approximant. We give a short 
introduction to all the considered post-Newtonian (pN) model waveforms. 

m ^ 

\ • The TaylorTl approximant [2] employs the instantaneous gravitational 

wave (GW) angular velocity w and £, the far-zone flux of the pN accu- 
rate energy E, for an evolution equation of both GW phase (j) and lj at 
time t. 
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dt (dS/dx) 

The TaylorT4 approximant (not used in LAL) is a slightly modified version 
of Tl and uses, as the only difference to Tl, a Taylor-expanded version of 
the ratio in Eq. ([3|). 

The TaylorF2 approximant is the Fourier domain pendant to TaylorTl, 
obtained by a stationary phase approximation |3I4) . 

The TaylorEt approximant [5] uses the orbital binding energy as the quan- 
tity to be evolved by radiation reaction. 
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where is the symmetric mass ratio rj = mi 7712/ (mi + 7712)^. 

For our studies, the post-Newtonian expansions of Eq. ([3]), GW luminosity £, and 
expressions in Eq. ((4]) are used up to 3.5pN order. AU the necessary references 
and exphcit formulas, partly only given with dots here, can be extracted from the 
investigation we like to summarize shortly [1] . 

2. Issues in GW phasing and data analysis 

For our data analysis considerations, the time domain GW forms are constructed 
in the following way. The orbital evolution is set to start at the initial instant tg 
when the instantaneous GW frequency of few kicks in the detectors bandwidth. It 
is necessary to transform few into the binding energy in a pN accurate way for the 
Et waveform and to let few or E evolve due to the model. The upper boundary for 
the GW frequency is that of the last stable circular orbit of a black hole having the 
total mass m of the binary system. The merger and ringdown phases are excluded 
from our analysis. 

3. Results and summary 

For the equal-mass case, the investigation showed up that for all three employed 
interferometer models, the systematic errors of the estimated masses continuously 
rise from to the order of magnitude of ^ 6 — 8% when going from m — 1.4 to 
m = 40 solar masses, supposing that the TaylorEt is the signal and TaylorTl and 
T4 are taken to be the templates. The fitting factors have a minimum at the mass 
in the signal of ^ IOA/q, which is always below 96%. Additionally, the FF values 
are always below 96% for advanced LIGO and Virgo for all the considered cases. 
For the unequal-mass case, the results behave qualitatively somewhat different. To 
give representative numbers, let us show the following table for 77 = 0.1875 (for the 
mi/m2 = 1/3 signal family) and advanced LIGO at 3.5pN order, presenting the 
FFs, the ratios of the estimated total mass with respect to the signal's one, and the 
estimated 77, as functions of template masses in the first row. 
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FF 


West /A/0 


Vest 


3:9 


12 


Tl 


0.98 


10.16 


0.250 






F2 


0.97 


10.23 


0.248 


5:15 


20 


Tl 


0.96 


16.91 


0.250 






F2 


0.97 


16.96 


0.250 


10:30 


40 


Tl 


0.98 


33.60 


0.250 






F2 


0.95 


33.82 


0.249 
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As it is clearly visible, our code evaluated the FF at large biases of the estimated 
total mass West and mass ratio rjcst ■ We cannot give a full overview here and like to 
refer to the more detailed presentation of the numerical results [T]. 

One result of this investigation is that fitting factors for > 95% for unequal- 
mass binaries can be only achieved at the cost of highly biased rricst and ?7ost • The 
faithfulness for the Tl and Et comparison turned out to be ^ 0.45. 

This investigation is only able to sketch the efficiency of capturing one model 
with the other. The question to be asked is what kind of analytic waveform is 
closer to reality can partly be answered by numerical relativity (NR) , which is only 
able to model few GW cycles at the end of the binary inspiral. It turned out that 
the TaylorT4 waveform is, for some cases, very close to NR. With respect to other 
waveforms that are closer to NR in the sense of matched filtering, it has unfavorable 
properties like artificial introduction of large parameter biases. This investigation 
confirms what has been presented in more detail in another recent work j6] and has 
been extended for TaylerEt to eccentric orbits [7|, where in the latter case we are 
left with the demanding question for numerical results for the eccentric inspiral. 
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